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UNCERTAINTIES IN DERIVED TEMPERATURE-HEIGHT PROFILES 

R. A. Minzner 

Atmospheric and Hydrospheric Application Division 
ABSTRACT 

Nomographs have been developed for relating uncertainty in temperature T (or 
in the corresponding scale height H) to uncertainty in the observed height profiles 
of both pressure p and density p. The relative uncertainty 5T/T is seen to de- 
pend not only upon the relative uncertainties 5 p/p or 5p/p, and to a small extent 
upon the value of T or H, but primarily upon the sampling -height increment Ah, 
the height increment between successive observations of p or p. For a fixed 
value of Sp/p, the value of 6T/T varies inversely with Ah. For Ah =\/2 -H, which 
is of the order of 10 km, 6T/T - 5p/p. For Ah = 0.01km, 5T/T is about 1000 
Sp/p, while for Ah = 35km, 5T./T is about 0.3 5p/p. In the case of T derived 
from density-height data, this inverse relationship between 5T/T and Ah applies 
only for large values of Ah, i.e. , for Ah > 35 km. For Ah < 1km, ST/T d 5 p/p, 
independent of the size of Ah. No limit exists in the fineness of usable height 
resolution of T which may be derived from densities, while a fine height resolu- 
tion in pres sure -height data leads to temperatures with unacceptably large 
uncertainties. 
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SYMBOLS 


Equation or 


Basic Quantities and Coefficients Figure 

A The dimensionless coefficient indicating the factor by 

which (Sp/p) is amplified to yield (5T/T) (20) 

G The standard geopotential gravitational constant 

9. 80665 m 2 sec~ 2 (m')" 1 (1) 

H Scale height in units of geopotential height, m' (or km' ) (1) 


h Geopotential height in units of m' (nearly equal to geo- 
metric height, but accounting for height variation of the 
acceleration of gravity: i.e., Gdh = g(z)dz (1) 

k A dimensional constant 0.0341632 K/m' resulting from 
the combining of the three constants GM/R into a single 


value (4) 

M mean molecular weight of air 28, 9644kg kmol" 1 (1) 

p atmospheric pressure Nm -2 (1) 

R Universal gas constant 8.31432 x 10 3 joules K' 1 kmol" 1 (1) 

S The value of a particular series of terms involving 

density ratios (30) 

T Atmospheric Temperature K (1) 

x A general designation for a function of a number of 

generalized variables (14) 


vii 



Equation or 

Basic Quantities and Coefficients (Cont’d.) Figure 

y A general designation for the set of variables of which 

x is a function (I 4 ) 

a A subscripted dimensionless coefficient with a value 



near unity 

(23) 

p 

atmospheric density kg m' 3 

(7) 

Single Subscripts 


1 

Refers to a particular height hj so that p* , T } , P* , 



and h } all refer to the same point in a. given equation 

(4), (H) 

2 

Similar to 1 

(4) 


Refers to an individual member of a set of re- 
lated variables having the general designation y ■ (14) 

Also refers to a variable in an isothermal atmosphere 
as P i in the graph of (p* /p r ) versus h in Figure 1 
A general designation for an integer which may vary be- 
tween 2 and q, and which is simultaneously associated with 
a geopotential height as h. and with the related density as p. (21) 

A general designation for an integer which may have any 
positive value and which is associated with the lowest 
density-height data point (i.e. , h q , p q ) involved in a parti- 
cular evaluation of an integral to determine the value of the 
related temperature T . (H) 

4 
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Single Subscripts (Cont'd.) 

q Also used as where q implies the number of terms 
in the series S 

r Designates a specific reference value for the basic 
quantity of which it is a subscript, i. e. , p r and p r 
Double Subscripts 

Refers to a particular height layer as between h 2 and h } , 
where hj and h 2 are ordered in such a way that the 
quantities within the subscripted parentheses have par- 
ticular signed significance, i. e. , Ah is always positive 
while AH is positive only when dT/dh is positive in the 
related layer 
Similar to 1, 2 
Similar to 1 , 2 
Similar to 1,2 

designating maximum value of S^, as in Figure 6. 
designating minimum value of S , as in Figure 6, 

H 

Operators and Functions 

Overbar as in T or H indicates the mean value of T or H 
for the related layer 
V Square root as V2 

/ Integral as /p(h)dh 


j-lj 

j»j+l 

q-l.Q 

q, max 
q, min 


Equation or 
Figure 


(30) 


Figure 1 


(4) 


(10a), (23) 
and (24) 

(25) 

(26) 
(27) 


(36) 


(4) 

( 20 ) 

( 11 ) 
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Operators and Functions (Cont’d.) 
d Partial differential as in dx 

d Differential as in dh 

£n Natural logarithm as £np and Cnp 
A An increment as in Ah and AH 

5 An increment or random uncertainty as in Sp [f 
5Ah, STj, 5 pj, 6T q , and 5p q 
2 A summation of a set of terms 


Equation or 
Figure 

(14) 

( 1 ) 

(1) and (7) 
(4) and (10a) 
(14) , (15) 
and (21) 
(12), (14) 
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AN ANALYSIS OF THE ERRORS ASSOCIATED WITH 
THE DETERMINATION OF ATMOSPHERIC TEMPERATURE 
FROM ATMOSPHERIC PRESSURE AND DENSITY DATA 

I. INTRODUCTION 

Relationships between pressure, temperature, and height in the earth’s atmos- 
phere are well known, and for many years have been the basis for height deter- 
mination in balloon- radio sonde flights making observations of pressure and tem- 
perature in meteorological probings of the lowest 30 kilometers of the earth’s 
atmosphere. When radar-tracked rocket vehicles extended the potential atmos- 
pheric probing capabilities to heights up to 100 kilometers, into regions where 
pressure is still measurable, but where existing technology does not allow for 
immersion sensing of temperature, these same mathematical relationships 
were used to extract temperature from the measured pressures and radar-derived 
rocket-height data. The unsmoothed temperature-height values from many of 
these soundings represented a very jagged height profile,, with the degree of jag- 
gedness apparently increasing as the height increment between successive data 
points decreased. Investigators usually have been unable to determine how much 
of this jaggedness represents real temperature variability, and how much is 
attributable to measurement error. 

The hydrostatic equation and the equation of state lead to another set of 
height relationships: those between atmospheric density, temperature, 
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and height, such that temperature may also be computed from density- 
height data. In cases where such computations have been made, partic- 
ularly in the height region of 30 to 100 km, the resulting temperature -height pro- 
files appeared to be less jagged than those derived from pressure-height data 
with comparable height resolution. 

The apparent difference in the jaggedness of density-derived temperature -height 
profiles from those associated with pressure-height data suggest that the height 
increments of the pressure and density data do affect the uncertainty in the de- 
rived temperatures, and that the influence of height increments in relation to 
density data may be different from that in relation to pressure data. Obviously, 
the error analyses presented in this paper, involving both the pressure- 
temperature -height relationship, and the density- temperature-height relation- 
ship, are needed. The error analyses confirm the fact that important differences 
do exist between these two sets of relationships, particularly in regard to the in- 
fluence of the height increment on the propagation of measurement uncertainties 
into the temperature-height profile. These differences strongly favor the use of 
density-height data over pres sure -height data. 
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II. FUNDAMENTAL CONSIDERATIONS 


A. Pressure-Height Relationships 

The equation of state, when combined with the differential form of the hydrostatic 
equation to eliminate density p, yields an expression frequently referred to as the 
hypsometric equation: 


d£np _ -GM _ -1 
~dh W H 


where 


p is atmospheric pressure , 

h is a measure of the height above sea level, in geopotential meters m' 
T is absolute temperature of the atmosphere at h, 

R is the universal gas constant, 

M is the mean molecular weight of air, 

G is a constant when h is expressed in geopotential, and 
H is the scale height in geopotential units. 


Solving equation (1) for T and H, respectively, yields 


_ GM -dh 

T — 

R dfinp 

and 


( 2 ) 


FT _ -dh 
GM dJ2np 


( 3 ) 


When values of finp versus h are known from numerical data rather than from an 
analytical function, it is convenient to replace the expression -dh/dBnp, which 
applies to a specific height, with a numerical approximation (h 2 -h T )/(.Cnp 1 -Cnp 2 ), 
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In this approximation, p 1 is the pressure at h- 1# and p 2 is the pressure at h 2 . The 
approximation therefore represents the mean value of the. reciprocal of the deriv- 
ative over the height interval (h 2 -hj). When the point value of the derivative is 
replaced by the numerical approximation, in both equations (2) and (3), the related 
values of T and H no longer apply to a single height, but rather become the mean 
values T and H, respectively, associated with the height interval (h 2 -hj). Thus, 


we have 


GM 


T = 


h 2 - hj 


= k 


(Ah) 1>: 


R £npj -Cnp 2 


Cnpj - £np 2 


( 4 ) 


where k = GM/R = .034163 °K/m’ , and (Ah) 1 2 = h 2 -h 1# and where, from the 
relationship between T and H implicity In equation (1), we see that 


H = 



(Ah) l 2 


( 5 ) 


2np 1 - £np 2 £np t - £np 2 

From equation (5), we obtain the following expression which will be of special 
importance in the uncertainty analysis in Section III: 


1 

Cnp 1 - gnp 2 



( 6 ) 


The interrelationships between Cnp, dCnp/dh, T, and H, all as a function of 


height as expressed by equations (1), (2) and (3), are shown picto rally by four 
of the six height- related graphs of Figure 1. This figure, which depicts the 
properties of a portion of the 1974 U. S. Standard Atmosphere (in preparation), 
was developed from a preliminary set of abbreviated tables representing a por- 
tion of this revised standard atmosphere (Kantor and Cole 1973). Geopotential 
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Figure 1. The U.S. Standard Atmosphere as Represented by Graphs of Five 
Atmospheric Properties versus Height, i.e. , Temperature, Scale Height, 
Natural Log of Density Ratio, Natural Log of Pressure Ratio, and Derivative of 
Natural Log of Pressure Ratio with Respect to Height, Plus a Sixth Property 
for Comparison, i.e. , Natural Log of Density Ratio versus Height for an 
Isothermal Atmosphere 
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d £ n p/dh (mV 1 TEMPERATURE (K) SCALE HEIGHT (km') 



height is scaled linearly along the abscissa of Figure 1. Atmospheric pressure; 
which is plotted in the form of the natural logarithm of the ratio <p/p r ), has its 
ordinate scale at the left of the figure, and is depicted by the lower of the three 
lines diagonally crossing the entire figure from upper left to lower right. In this 
presentation of pressure, p r in the ratio (p/p r ) has been made equal to one Newton 
per square meter (i.e. , 1 N m -2 ), so that the numerical values of (p/p r ) are 
those of p inN m' 2 , while the scale of Cn(p/p r ) at the left of the figure implies 
no dimensions. 

The graph of d£np/dh, which quantity is negative and inversely proportional to T, 
shows that the slope of Cn(p/p r ) versus h, or equivalently the slope of Cnp versus 
h, is constant only in height regions over which T is invariant. These portions of 
the graph of dfinp/dh versus height corresponding to constant T (i.e, , the three 
height intervals 11 to 20 km' , 47 to 57 km' , and 84. 852 to 89. 716 km/ ) consist 
of horizontal straight-line segments, while the remainder of that graph consists 
of 5 curved segments. The scale height, which is directly proportional to T, is 
seen to remain within the range of 5,45 and 7. 95 km' , with a median value of 
6.5km'. < 

B. Density-Height Relationship 1 

When the differential form of the hydrostatic equation is combined with the differ- 
ential form of the equation of state to eliminate the differential of pressure, we 
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obtain the following relationship between temperature and the density- height 
function: 


-dh 
d £ np 



( 7 ) 


In the absence of an analytical expression for p(h), a numerical approximation 
for dh/dKnp leads to a mean value of temperature T for the height interval h 2 - hj 

( 8 ) 


T = 


h 2 h 1 


Cnpj - £np 2 


GM 


T a- T , 


R 


h 2 - h l 




In terms of scale height, this expression becomes 


H = 


h 2 -hi 


2npj - 2np 2 


1 + 


H 2 -Hi 
— 1 

h 2 ' h l 


(9) 


It is important to note that while equation (8) expresses the mean temperature 
for the layer h x to h 2 > the expression involves the average gradient of T with 
respect to h within that layer. A similar situation prevails for the expression 
of the mean scale height in equation (9), where H = (H 2 + H^/2. Since prior 
knowledge of the value of neither (T 2 - T 1 )/(h 2 - hj) nor (H 2 - Hj)/(h 2 - hj) is 

i 

generally available, equations (8) and (9) are not of themselves useful relation- 
ships for temperature-height determination. One version of the following ex- 
pression derived from equation (9), however, will be of special importance in an 
uncertainty analysis in Section IV. This expression is 

1 (H-+HJ/2 

— = — GO) 

Cnpj - £n p 2 (h 2 - hj ) + (H 2 - Hj ) 

Equation (10) is written for the case in which h 2 > h lt as when data are being 
analyzed from lower to greater heights, such that p 1 > p 2 . If the equation is to 
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apply to data being analyzed from greater to lower heights, as is the preferred 
case for the integral form of the density- height function, involving the normal 
atmosphere, h 1 is greater than h 2 such that p 2 > p 1 , and equation (10) must be 
rewritten as 

1 O^+Hp/2 / H 

Cnp ? - Cnpj (hj -h 2 ) + (Hj ~H 2 ) \Ah + AH 

In this expression Ah is always positive, and AH is positive in regions of positive 
temperature gradients, zero in regions of zero temperature gradients, and neg- 
ative in regions of negative temperature gradients. 



C. Density- Height Relationship 2 

The integral of the hydrostatic equation, when combined with the equation of state 


to eliminate the pressure, yields the following directly useful integral expression 
relating temperature to density: 


Pn R P n A 


p(h) dh 


(ID 


'q “ q 

The temperature T q is not the mean value of temperature for some layer, as in 
the case of f in equation (4), rather, it is the temperature at the specific height 


It has been demonstrated (Minzner, etal. 1964 and 1966) that the properties of 
equation (11), when applied to a helium atmosphere, differ markedly from those 
observed when the equation is applied to an argon atmosphere. With an air at- 
mosphere having a mean molecular weight of about 29, the properties of this 
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equation are similar to those found when the equation is applied to an argon at- 
mosphere in which the integration optimumly proceeds from the greatest to the 
lowest altitude of the density-height data. In this situation, T t , p t , and h l are 
each associated with the greatest height of the data set, while T , p , and h a are 
each associated with the running value of h. This is the height for which the value 
of T q is being computed, a height which varies progressively from h L to the low- 
est height of the data set as the calculation of the profile proceeds. 


Because p(h) is usually known numerically rather than analytically, it is conven- 
ient to replace the integral of equation (11) with an appropriate series approxima- 
tion, one of which is governed by the logarithmic trapezoidal rule (Smith 1965, 
Minzner, etal. 1965). With the use of this approximation, equation (11) may be 
rewritten as 


T 


q 



T 1 + 


GM 

~R 


< y* (hj _ 2 -hj) (pj - pj . 
Oq J = 2 ftnpj - Snpj _ x 


( 12 ) 


The logarithmic-trapezoidal-rule approximation is a particularly suitable one 
for the integration of atmospheric density with respect to height, because, in a 
graph of measured values of finp versus h, the straight-line segments between 
successive data points exactly represent the conditions of the logarithmic trape- 
zoidal rule, and also very closely represent the conditions of the real atmosphere. 
The latter is evident from the nearly straight-line segments of an idealized ver- 
sion of the real atmosphere shown in the graph of £n(p/p r ) versus height in Fig- 
♦ ure 1. For this graph, which represents the densities of the 1974 U.S. Standard 
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Atmosphere, the value of p r was chosen to be 1 x 10“ 5 kg m -3 , so that the single 
ordinate scale at the left of the figure applies to both fin(p/p T ) and ftn (p/p r ). The 
shape of the curve £n(p/p r ) versus h is identical to the shape which the curve Cnp 
versus h would have, since the values of £n(p/p r ) vs. h and of 2np vs. h are off- 
set by the constant difference £ np r at all heights . , 

Figure 1 also contains a graph of £n(p t / p T ) in the form of a continuous straight 
line immediately above the graph of Cn(p/p r ). This single straight line 
diagonally across the entire figure is characteristic of the density of an isother- 
mal atmosphere extending upward from a height of 7. 625km where the standard- 
atmosphere density is 5.8281 x 10" 1 kg m' 3 , and where the temperature and the 
corresponding scale height, are 238. 587 K and 7, 000 km' , respectively. For this 
specialized atmosphere, the logarithmic-trapezoidal rule represents a series 
which exactly duplicates the integral of equation (11). The small deviations of 
the slopes of £n(p/p r ) versus height from those of %n(p { /p r ) versus height show 
the small influence of the variation of atmospheric temperature from the fixed 
value, 238. 587K, (in the height region of 10 to 90 km' ) on the general shape of 
the curves of £n(p/p r ) and £np versus h. (The influence of temperature variation 
upon the curves of £n(p/p r ) and £np versus h is similarly small. ) 

In the evaluation of T by equation (11), or by any appropriate approximation of 

4 

that equation as exemplified by equation (12), a knowledge of the initial tempera- 
ture Tj is of importance only for the upper regions of the profile. For h q = h it 
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the value of (p 1 /p q )T 1 is exactly T { , since p { = p q and the integral term is, of 
course, zero. As the value of h q decreases from h lf i.e. , as(hj - h q ) increases, 
the value of the relative contribution of the density- ratio term decreases, while 
that of the integral term correspondingly builds up. When (h L - h q ) increases 
from zero, first by one scale height and then by three scale heights, the relative 
contribution of the density-ratio term to T q decreases from 100% to about 37%, 
and then to about 3%, because of the large decrease in the value of (/^/pq) over 
these height regions. Simultaneously, the value of the integral term grows cor- 
respondingly. At heights of more than three scale heights below h lf the value of 
T q is determined almost completely by the integral term alone. Figure 2 shows 
the relative contributions made to T by the density-ratio term, and by the inte- 
gral term of equation (11) as a function of the range of the limits of integration 
(h 1 - h q ) expressed in units of scale height. 

The graphs in Figure 2 are based on the density-height profile for an isothermal 
atmosphere (T = 238. 587 K) for which £np is a linear function of height as pre- 
viously shown in the graph of Cn(p i /p r ) in Figure 1. Only small variations from 
the values of the two terms of equation (12) depicted in Figure 2 would be seen 
for calculations based on a real, variable -temperature atmosphere. 


The elimination of T between equation (12) and equation (3) yields the following 
expression for scale height: 


H 


q 




A (hj - t ~ hj) (Pj - Pj . t ) 


(13) 
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RELATIVE CONTRIBUTION TO T q 



Figure 2. Relative Contribution Made to the Calculated Temperature by the 
Density-Ratio Term and by the Integral Term of Equation (IX) as a Function of 
the Size of the Altitude Range of Integration, for a 238 K Isothermal 
Atmosphere 

The relative contribution of each of the two terms of this equation to H q , as a 
function of the range of integration (hj - h q ), follows exactly the same pattern 
shown in Figure 2 with regard to T q . 

Atmospheric soundings of both pressure and density have for many years served 
as the basis for the determination of the height profiles of temperature and scale 
height. From pressures, these profiles are obtained through equations (4) and 
(5) respectively, while from densities, they are obtained through one or another 
version of equations (12) and (13), respectively. It is obvious that the uncertainties 
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of these derived temperatures and scale heights are a function of the uncertainty 
in the measured quantities, pressure or density. Somewhat less obvious is the 
fact that the height interval between successive observations of pressure or 
density strongly influences the propagation of the observational uncertainties 
into the computed temperatures and scale heights, A rigorous error analysis 
of both of these pairs of equations demonstrates this situation. The pressure- 
related equations are analyzed in Section in, while the density- related equations 
are analyzed in Section IV. 
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III. UNCERTAINTY of temperatures deduced from pressures 

The error-analysis method employed is the first -order Gaussian method, wherein 
each variable y * entering into the expression of a particular function of these 
variables x(y t ) is assumed to have an observational uncertainty fiyj which meets 
the conditions of a Gaussian or normal distribution about y. , where y i is the 
mean of a set of individual observations of the ith variable, or the true value 
of the ith variable. Thus, if the value of x is determined from the functional 
expression x(y. ), the value of fix, the implicit uncertainty in x, is given by 

(6x) 2 =£ • 8y^ (14) 

Applying this relationship to the variables of equation (4) yields 



It may be shown that the partial derivative of T with respect to any one of the 
independent variables in equation (4) has the general form 

df - kAh 5y: 

S yi ( 16 ) 

dy { (finpj-finpp 2 yj 
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where Ah is associated with (h 2 - hj) without benefit of subscripts. Applying 

this generalized expression to equation (15) yields 
_ +kAh 


ST = 


(finpj - £np 2 ) 2 j_\ Pj 


\ 2 /SpA 2 /5Ah\ 
— ' + — + — - 
\ p 2 ) V Ah / J 


'A 


( 17 ) 


where 5p 1 and Sp 2 are the pressure uncertainties of two consecutive pressure- 
height values, and 6 Ah is the uncertainty of the height interval Ah between 
the corresponding two pressure-height values. 

Dividing each side of equation (17) by the appropriate side of equation (4) yields 
the relative uncertainty 


ST 


1 


-p Cnpj - £np 2 


'*p. N 




5 Ah' 

Ah j 




( 18 ) 


Pi/ \P2/ 

It is convenient to assume that all the uncertainty of a pres sure -height point is 
in the pressure, and that the height increment is exact. In this case, the values 
of both Sp L and Sp 2 are correspondingly increased over their actual values, and 
the term involving SAh in equation 18 vanishes. It is further assumed that the 
relative uncertainty in measured values of p is the same at hj and h 2 , i. e. , 
^Pi/Pi - ^p 2 /p 2 = 5 p/p. These two assumptions permit equation (18) to be re- 
duced to the simpler form, 


ST 


V? 


8p 


T - Cnp 2 p 

From equation (6) it is apparent that 1 /(£np { - £np 2 ) in equation (19) may be 
replaced by (H/Ah)j 2 where H is the mean scale height for the layer 


( 19 ) 
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k 2 - h 1 = (Ah)j 2 . Thus, equation (19) may be rewritten as 



This equation indicates that (5T/T), the relative uncertainty in T the mean 
temperature of a layer (Ah) 1 3 (i. e. , the layer bounded by pressures p, and p 2 ), 
is equal to the relative uncertainty in the pressure measurements (§p/p) times 
a coefficient A i 2 . This coefficient is seen to depend upon only two quantities, 
the mean temperature T of the layer h 2 - h 1 (T being implicit in H), and (l/Ah), 
the reciprocal of the thickness of the pressure -sampling interval (h 2 - h L ). Since 
the maximum value of T in the earth's atmosphere below h - 100 km' is less 
than twice its minimum value, and since Ah can, in principle, be made to vary 
over many orders of magnitude, it is apparent that Ah is the dominant factor in 
determining the propagation of (Sp/p) into (5T/T). 

A graph of the coefficient A } 2 as a function of Ah for three specific values of 
T, 169.10 K, 241.57 K, and 314.04 K, is given in Figure 3. The first and 
third of these temperatures were selected in part because they are close to the 
lowest and highest atmospheric temperatures normally observed at heights below 
120 km, while the second is the mean of the extremes. In addition, these three 
values were specifically selected to correspond, respectively, to a particular set of 
three scale heights,' i. e. , (7/V2) , (lO/v/2), and (13/V2) km' . In each of these three 
cases, the value of A 12 is unity when Ah (equal to \/2"* H) is a particular integer 
multiple of one geopotential kilometer, i.e., 7, 10, and 13km', respectively. 
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VALUE OF COEFFICIENT ^ (dimensionless) 



T = 31 4.04 K 


T= 169.10 K 


UNITY 



PRESSURE-SAMPLING HEIGHT INTERVAL Ah (m‘ 


Figure 3. Value of the Coefficient ,f A { ” in Equation (20) as a Function of 
Pressure -Sampling Height Interval for Each of Three Atmospheric Tempera 
tures Representing Approximately the Maximum, Median, and Minimum 
Values, Respectively, for Heights Below 100 Geopotential Kilometers 





It is interesting to note that, for the entire range of normally observed temper a- 
tures at heights below 100 km \ the pressure-sampling interval Ah corresponding 
to unity for the coefficient A 12 varies between the limited range of 7 to 13km'. 
The pressure -sampling height interval, however, may actually vary over 
several orders of magnitude depending upon the design of the measuring system. 
The value of the coefficient A x 2 could correspondingly vary over several orders 
of magnitude depending upon the choice of Ah. 

Concentrating on the median value of T, and allowing the pressure -sampling 
interval Ah to decrease, first fromlO 4 m' tolO 3 to!, and then tolO 2 m , causes the 
value of A, ? to increase first from ItolO, and then to 100, respectively, such 
that the uncertainty in the value of T as expressed by equation (20) increases by 
identical factors for a fixed uncertainty in the pressure. Conversely, as Ah is 
increased from 10 4 m'to 10 5 m', the value of A 1>2 decreases from unity to 0. 1, 
such that the uncertainty in a related mean value of T decreases correspondingly. 

Figure 3 shows that when mean temperatures are deduced for successive layers 
from pressure observations having a fixed uncertainty, finer height resolution 
in the resulting temperature -height profile is obtained at the expense of tempera- 
ture accuracy, while increased temperature accuracy is achieved at the expense 
of height resolution. Since the ultimate accuracy of pressures obtained from 
any high-altitude balloon-borne or rocket-borne pressure -sensing device is 
limited by considerations of various perturbing phenomena such as outgassing, 
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boundary layer, and shock wave, the value of 5 p/p has a practical lower bound. 
Thus, ( 6T/T) • (AhA/IT • H) which is equal to 5 p/p similarly has a practical 
lower bound, and for this minimum value of 5 p/p, 6T/T is governed by the 
value of the ratio (Ah/H), 

Figure 3 shows the effect of variations of both T and Ah on the value of the 
coefficient A : 2 = {\/2 * H/Ah). By expressing Ah in multiples of one scale 
height, die number of variables in equation (20) is effectively reduced by one, 
and values of 5T/T can be plotted as a function of Ah = n ♦ H (where n is the 
value along the abscissa) for any particular value of 6 p/p. Figure 4 presents 
such a graph for each of nine values of pressure uncertainty. 

In Figure 4, the straight-line graph for 1% uncertainty in p intersects the co- 
ordinate for 1% uncertainty in T at the abscissa coordinate value of 1.414. 
This same straight-line graph for 1% uncertainty in p intersects the ordi- 
nate values of 10%, 100%, and 1000% at abscissa values of Ah equal to 
0.1414, 0.01414, and 0.001414, respectively. This series of four suc- 
cessively decreasing values of Ah, each being one tenth of the preceeding 
one, yields a geometric series of four increasing values of 5T/T, each being 
ten times the preceeding one. A similar situation prevails for each of the eight 
other values of 5 p/p, for which lines have been plotted. The value of 6T/T is 
obviously varied over orders of magnitude by comparable variations of the value 
of the ratio H/Ah. For height increments with a value of>/Z times one scale 
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UNCERTAINTY IN MEAN TEMPERATURE (percent) 



PRESSURE-SAMPLING HEIGHT INTERVAL Ah (units of one scale height) 

F ignre 4 . Percent uncertainly in the Mean Temperature of an Atmospheric Layer 
(for any one of nine percent uncertainties in the pres sure -height data) as a func- 
tion of the Pressure -Sampling Height Interval (i.e. , the thickness of that Layer) 
when that Mean Temperature is Deduced from the Atmospheric Pressure at the 
Upper and Lower Boundaries of that Layer » and When the Mean Temperature is 
At or Near 241. 57 K. 



height, relative temperature uncertainties are seen to be equal to relative pres- 
sure uncertainties for all nine values plotted. 

In order to show simultaneously the small additional influence produced by the 
allowable range in the value of T as it appears both in 5T/T as well as intrinsi- 
cally in H, the single-line graph for each value of 5 p/p in Figure 4 is expanded 
into a band in Figure 5 where the pressure -sampling height is expressed in 
meters as in Figure 3. 

The lower left-hand edge of each band corresponds to T - 169. 10 K, while the 
upper right-hand edge of each band corresponds to T - 314. 04 K with other 
points across each band corresponding linearly to intermediate temperatures. 
Thus, Figure 5 serves as a nomograph from which one may estimate the per- 
cent uncertainty in particular atmospheric temperatures computed from pressure- 
height data measured with a specified uncertainty over particular height 
increments. 

Three sample applications of this nomograph are cited: In the first it is desired 
to determine the maximum pressure-gauge uncertainty allowable to achieve a 
1% uncertainty in a mean temperature of about 240 K for a layer thickness 
of 100m'. We look for the intersection of the 100m' abscissa value with 
the 1% ordinate value, and find that it lies near the 242 K value of the 0.01% 
pres sure -uncertainty band. Thus, pressures would have to be measured with 
an uncertainty no greater than 0. 01% to achieve the desired results. Since it is 
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UNCERTAINTY IN MEAN TEMPERATURE (PERCENT) 


PRESSURE-SAMPLING HEIGHT INTERVAL Ah (m 


Figure 5. Nomograph for Estimating Relative Uncertainty in Temperatures 
(Calculated from Pres sure -Height Data) as a Function of the Pressure -Height 
Sampling Interval, and as a Function of the Value of the Temperature, for Each 
of Nine Values of Relative Uncertainty in the Pressure -Height Observations 



essentially impossible to achieve such a small uncertainty in any rocket or bal- 
loon measurements of pressure, this combination of height resolution and 
temperature uncertainty is essentially impossible to achieve from pressure- 
height data. 

In the second example, we assume a pressure-gauge uncertainty of 3 percent, 
and a mean temperature of 180 K, and seek the temperature uncertainty associ- 
ated with particular sampling-height increments. If the pressure sampling- 
height interval is 3 kilometers, this ordinate value is seen to intersect the 
appropriate region of the band for 3 percent pressure uncertainty at a value 
corresponding to a temperature uncertainty on the ordinate scale of about 7.4%, 
or about 13. 3 K. A doubling of the layer thickness to 6 km would halve the 
temperature uncertainty. 

The third example involves the inverse problem of estimating the uncertainty 
in the computed pressure differential associated with an assumed isothermal 
layer of fixed thickness. If the layer has a thickness of 3 km, as may be the 
situation in the grenade experiment (Nordberg and Smith 1964), and if there is 
an uncertainty of about 2% or 5 K in an assumed mean temperature of 250 K, 
the graph shows that the combined uncertainty of the two boundary pressures 
lies between 0.3% and 1. 0% or about 0. 6% on a logarithmic scale. For a layer 
thickness of lkm, the boundary-pressure uncertainties would increase by a 
factor of 3. If the assumed isothermal layers were used to generate 
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density-height profiles rather than pressure-height profiles, the uncertainties 
in the densities would be somewhat larger than in the pressures because of tlie 
required vector addition of the temperature uncertainty and the pressure un- 
certainty to obtain the density uncertainty. 
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IV. UNCERTAINTY IN TEMPERATURES DEDUCED FROM DENSITIES 
Applying the first-order Gaussian method, equation (14), to the determination 
of the uncertainty in temperatures deduced from density -height data through 
equation (11) leads to 



provided that we assume all the uncertainty in a density -height point to be con- 
centrated in the density. Because scale height combines the temperature with 
several constants, its use facilitates the expansion and analysis of an equation 
like equation (21). Consequently, it is convenient to rewrite that equation as 



Both equations (21) and (22) involve (q + 1) terms, where q is the number of 
density-height data points. The first term is associated with the uncertainty 


of the temperature or the scale height at h lt while each of the q additional terms 
deals with the uncertainty of one of the successive q density-height data points, 
respectively. Three of the (q + 1) terms have a unique format. These are the 
first, the second, and the last terms; those involving ST 1 (or SIty), 5p 1 , 
and 5p q , respectively. The remaining terms, however, the third through the 
qth term involving bp 2 through 5p q .j , respectively, have a common format. 

Thus the sum of these comm on -format te nis may be expressed as the summa- 
tion of a general term. Applying the operators indicated in equation (22), dividing 
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both sides of the resulting equation by H , and introducing equation (10a) yields 

* H k 

* 

the following equation: 



where 


a 


/H Ah \ 
ll2 ~ U ' Ah + AhJj 2 


/ H Ah \ 
VllJ = Ah + AH^.^J 


°ii + i 


/ H Ah \ 

\H q ‘ Ah + AHj j)j + 1 

/H Ah \ 

Vl ’ q = 


(24) 

(25) 

(26) 
(27) 


Equation (23), which is developed in Appendix A, is very much more complicated 
than the comparable uncertainty expression in terms of pres sure -height data as 
represented by equation (18). The considerable difference in the complexity 
between the two equations stems largely from the fact that equation (18) involves 
the uncertainty of only two pressure-height data points and the mean temperature 
for the corresponding height layer, while equation (23) involves the uncertainty 
of many consecutive density-height data points, each of which contributes at 
least a small amount to the value and uncertainty of the temperature or scale 
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height at height h q . It is emphasized that the quantities T q or H q in equations 
(11) through (13), and again in equations (21) through (23), represent values for 
specific heights, and not mean values for specific layers as in the case of the 
pressure-related equations. 

Each element of equations (23) through (27) is associated either with a specific 
height or with a specific layer. Elements having a single subscript, e.g. , H 1S 
Pi * Pj * Pj . i » Pj + i » Pq - 1 * P q ? and Hq signify the value of the particular quantity 
at heights h 1? hj, hj , v hj +I , h q _j, and h q , respectively. Quantities with a 
double subscript, e.g. , or 1 2 , l j, etc. , represent a quantity associated with 
particular layers, i.e. , the layers bounded by h 1 and h 2 , or by hj _ 1 and hj, etc. , 
respectively. Thus, n J 2 » as expressed by equation (24), represents an algebraic 
expression of four different quantities three of which, namely H, Ah, and AH, 
are associated with the layer h, to h 2 . In equation (24), H is the mean scale 
height for the layer Ah = hj - h 2 , while AH represents the change in scale 
height within that layer, i. e. , H* - H 2 , such that only in a region where the 
gradient of H or T with respect to height is positive will AH be positive. The 
quantity H q is, of course, the value of H at height h q . Similarly, o^ t . , a . 9 

and a q-i, q » as defined by equations (25), (26), and (27), respectively, each 
represent quantities associated with the particular appropriate layers* Equation 
(23) also includes the doubly subscripted quantities 

I . 

j. j + i 


H 

Ah + AH 


Ah + AH/j 2 ’ \Ah + AH/j „ \ ^ \Ah + AH 
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each of which represents one factor' of the right-hand side of equations (24), 

(25), (26), and (27), respectively. The influence of non-zero gradients in thfe 
temperature-height profile is impressed upon 8H q /H q through all the double 
subscripted quantities in equation (23). 

In order to more clearly see the influence of sampling-height interval alone on 
SH n , it is convenient to assume an isothermal atmosphere, and to correspondingly 

T. 

simplify equation (23). In this case, scale height would not vary, but would 
remain fixed at a value H over the entire region of integration, and the following 
relationships would apply: Hj=H q = H, H = H, and AH is zero In all layers. 
Under these conditions, equations (24), (25), (26), and (27) become unity, as 
does the ratio Hj/Hq, and equation (23) may be rewritten as 



It can be shown that 



28 



such that these two terms in the coefficient of (dp Q /p a ) in equation (28) cancel 
each other. Then, if we impose the additional condition that the height incre- 
ments between successive density-height data points are constant over the entire 
region of integration, this condition plus equation (29) permit the further simpli- 
fication of equation (28) to 



If we impose still another restriction, i.e. , that the relative uncertainty of the 
density data has the constant value dp/p for all heights within the range of inte- 
gration, 8p 1 /p 1 = SPj/pj = 5pq/ pq = 5p/p, and we may rewrite equation (30) as 



The restriction permitting this simplification, while somewhat unrealistic from 
the point of view of any measuring system, is acceptable because, it will be shown 
below, the relative uncertainty of only p q enters significantly into the value of 
5Hq/Hq. It is convenient to represent the sum of the series of the squares of the sev- 
eral density ratios in the coefficient of ( Sp/p ) 2 in equation (31) by (S ) 2 , such that 

H 



with this simplification equation (31) may be rewritten as 
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P i 
Pn 
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(33) 
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Finally, recalling from Figure 2 that the ratio ( p x /p q ) causes the contribution 
of H : to H q in equation (13) to become negligible when the integration has pro- 
ceeded downward from hj by more than about three scale heights, we recognize 
that the contribution of 6H 1 /H 1 to 8H q /H q in equation (33) must similarly be 
negligible for h q sufficiently below h 1 , Thus, for this condition, equation (33) 
may be approximated by 



(34) 


Also, since 5H q /H q is identically equal to 8T q /T q , it is convenient at this point 
to return from scale-height notation to temperature notation, and equation (34) 
is rewritten as 



(35) 


This equation represents the relative uncertainty in T q , the temperature at 
height h q , as deduced from a set of density-height data measured with a constant 
relative uncertainty (5p/p) over the entire range of a height region h 1 down to h q , 
where (h l - h q ) is equal to at least three scale heights, in a portion of a planetary 
atmosphere which is isothermal, at some value associated with H. These con- 
ditions might apply' to the earth’s atmosphere at heights above about 400km. 

Since equation (35) excludes several of the variables included in the most general 
form of the uncertainty expression as given by equation (23), this simpler version 
is more readily examined for the influence which variations in the sampling-height 
interval Ah alone have upon 5T q /T q . 
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Equation (35) for ST /T , in terms of the relative uncertainty in density, is very simi- 

4 4 

lar to equation (20) for ST/T in terms of the relative uncertainty in pressure, except 
that Sq in equation (35) takes the place of \^2 in equation (20). Therein lies a great dif - 
ference, because the quantity S q is a dimensionless function of three variables, H, 

Ah, and(h 1 -h q ), and has a value which varies between VJ and zero, asAhvai’ies 
from large values to zero. For Ah > 5 H , the value of S q approaches \/2asymptotically 
as Ah approaches infinity , independently of the value of H or (h L - h q ) . For Ah < 0. 2 H, 
the value of S q decreases linearly with decreasing values of Ah, for any fixed value of 

H or(hi -hq). 

The influence of variations of both Ah and (h 1 - h q ) upon S q , and implicitly upon 
STq/Tq, is seen in Figure 6. Here, for a fixed value of H = (10/>/2) km, a band 
of values of S q is plotted as a function of Ah on a fully logarithmic scale, with 
ordinate values of S q given on the right-hand side of the figure. The upper limit 
of the band is designated S q max , and is plotted as a line of long dashes, while 
the lower limit of the band is designated S qmin , and is plotted as a solid line. 

The range of the values of S q within the band, at any particular value of Ah, rep- 
resents the influence of variations of (h l - h q ) at that value of Ah, with the value 
of S q increasing as (h t - h q ) decreases. The smallest possible value of (hj - h q ), 
at any particular value of Ah, is of course Ah, and hence S q max depicts the locus 
of the values of S q associated with the condition that (hj - hq) = Ah. Thus, this 
upper limit of S q represents the value of equation (32) when that equation involves 
only the first two density-altitude points of a data set, i,e. , the value of S for 

4 

the case when 5T q /T q is being computed for the height h q associated with the 
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Figure 6. Graphs of Three Quantities as a Function of Density- Sampling Height 
Interval Ah, i. e. , (1) the Limited Band of Values of a Particular Series of 
Terms Defined by S q , and Related to Ah, (2) the Value of the Last Term of that 
Series, and (3) the Value of the Ratio of Mean Scale Height to Ah 


second-highest density- height data point in the sounding. For this situation, the 
right-hand side of equation (32) reduces to [\/2(p q - p, )/p q j . As the height for 
which one computes ST q /T q is lowered from hj , such that (h L - h q ) becomes 
increasingly large, the value of S q slowly decreases, and asymptotically ap- 
proaches a lower limit S q min when (h, - h q ) is greater than three scale heights. 


It is convenient to define the ratio of S q max to S qmin as R, and to examine the 
variation of R with respect to variations in Ah. As long as Ah is smaller than 
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about 0. 2 H, R retains a nearly constant value (a value which is seen graphically 
to be approximately equal to \/2), and both S qrnax and S qrnin are seen to make 
an angle of about +45° with respect to the abscissa. Because of the feature of 
a nearly constant slope of +1 for both Sq imax and Sq jm j n in the region where 
Ah < 0. 2 H, this region is hereafter designated as the small-increment regime. 

As Ah exceeds 0.2H, and increases toward 5H, the slopes of both S q max and 
S q min gradually decrease toward zero. However, the slope of S qmax decreases 
more rapidly than that of S q m ^ n so that S q max approaches S q min at a common 
value \/2, while the value of the ratio R decreases from V2 toward unity. For 
Ah > 5H, the value of R is essentially constant at unity, and the slope of the line 
common to both S qmax and S qmin is essentially zero. The feature of an essen- 
tially zero slope for this common line representing S q versus Ah, specifies a 
region which, because of the related large value of Ah (i.e. , Ah > 5 H) is here- 
after referred to as the large-increment regime. The region between the small- 
increment regime and the large-increment regime, i. e. , the region for which 

0. 2 H < Ah < 5 H, is hereafter referred to as the transition regime. 

/ 

In the large-increment regime, where S q niax and S q min are essentially identical, 
it is apparent that S q is essentially independent of (h t - h q ). In the small -increment 
regime, the concern for the influence of (h t - h q ) upon S q remains only for values 
of (hj - h q ) less than about 3 H. Because equation (35) is based upon the restric- 
tion that (h l - h q ) is greater than 3 H, and because S q approaches S q>m j n for 
such values of (h 3 - h q ) in all three regimes, it is immediately apparent that the 



general factor S q in equation (35) should be replaced by the specific value S q;fein . 


It is somewhat less apparent that 8p/p, the general expression for relative un- 


certainty in equation (35), aught to be replaced by $P q /P q , the specific relative 
uncertainty for h q , particularly when equation (35) is applied to the small- 
increment regime. The reasons for the latter replacement stem from the facts 
that the following dual relationship exists for the dual conditions (h x - h q ) > 3 H, 
and Ah < 0. 2 H: 

(1) The value of the series S q , which becomes essentially S q>m j n for 

(h l - h q ) > 3 H, simultaneously approaches the value of (P q - Pq-i)/P q > 
the last term in equation (32), defining S q . This near equality between 
s q min and (Pq “ Pq-i )/P q in the small -increment regime is demonstrated 
graphically in Figure 6. 

(2) The term (P q - p q-1 )/p q is properly associated with the specific un- 
certainty 8p q /p q rather than with 5p/p, This situation is evident from 
the fact that the quantity (p q - Pq-i )/p q I s a coefficient of Sp q /p q in 
equation (30), which, except for an isothermal condition, is a general 
expression equally applicable to both SH q /H q and ST q /T q . Thus, at 
least for the small-increment regime, equation (35) may be rewritten 
as 



( 36 ) 


Even in the large -increment regime, where the value of S qmin approaches 
[V2(P q - P q -i)/Pq] , the value of S qqTlin is still dominated by (p q - p q ^ )/p q . 
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(which is equal to unity in this regime), and therefore, S q min should still be 
associated with 5p q /p q . 

Because (p q - )/p q approaches unity in the large-increment regime, while 
Sq min approaches \f2, it is reasonable to write the following special form of 
equation (36) for that regime: 



This equation is seen to be analogous to equation (20) involving uncertainties in 
pressure-height data. 

In addition to the graphs of S q>max , S qmin , and (p q - P q _ } )/p q versus Ah, 
Figure 6 also contains a graph of the function (H/ Ah) versus Ah, for H = (10/ 
V%) km, or T = 241.57 K. This function, which obviously varies inversely with 
Ah, is seen to have a constant negative 45° slope for a fixed value of H. When 
this function is multiplied by S q min , a quantity which according to Figure 6 
is seen to have a +45° slope in the small-increment regime, the product for 
the small-increment regime is essentially a constant (independent of Ah) with 
a value near unity, for all values of H or T. This situation is depicted in Figure 
7 where the value of the product [S qmin • (H/Ah)l versus Ah is plotted in all 
three regimes for each of three values of H, (7A/2), (10/%/2), and (13A/2) km’ , 
consistent with the mean temperatures 169.10, 241.57, and 314.04 K, respec- 
tively. Because of the characteristics of this product, equation (35) as applied 
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Figure 7. Value of the Product of Two Factors (Comprising the Coefficient of 
the Uncertainty of Observed Densities) in Equation (36) as a Function of Density- 
Sampling Heig'ht Interval, for Each of Three Values of Mean Scale Height Con- 
sistent with Three Specified Temperatures 


to the small-increment regime may be replaced by 


ST n bo 


(38) 


Thus, three different expressions represent the value of 5T^/T q as a function 
of uncertainty in density-height data, equation (38) in the small-increment 
regime, equation (37) in the large-increment regime, and equation (36) in the 
transition regime. 


The application of the data from Figure 7, as representative of the three equa- 
tions noted, to each of nine particular values of uncertainty in density-height 
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data, yielded the nomograph shown in Figure 8. This figure shows the percent 
uncertainty in the temperature, derived from density -height data through equa- 
tion (12), as a function of density- sampling height interval for a band of normal 
atmospheric temperatures, 169 K to 314 K, and for a wide range of density- 
height uncertainty, 0.01% to 100%. This nomograph can serve as the basis for 
estimating the relationship between the uncertainty in each value of a set of 
density -height data to each point of the related Temperature-height profile, or 
vice versa, for the region of the earth T s atmosphere below about 120 km. 

The examples of hypothetical use of the pressure-data nomograph, discussed 
in Section III, could be readily transposed to apply to this density-data nomo- 
graph. In particular, Figure 8 shows that, in order to have no more than a 
1% uncertainty in a 245 K temperature, for a temperature height-profile with 
a height resolution of 100 meters, the relative uncertainty in the density obser- 
vations need not be smaller than 1%. This value is 100 times larger than the 
0.01% uncertainty required of pressure data to meet the same conditions, and 
represents a very significant relaxation of the measurement requirements. 

The estimates deduced from the use of this nomograph are theoretically accurate, 
as long as the data to which the nomograph is applied comply with the three re- 
strictions imposed during the preceeding development. It will be seen, however, 
that one of these restrictions may be essentially eliminated, while a second has 
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Figure 8. Nomograph for Estimating Relative Uncertainty in Temperatures 
(Calculated from Density-Height Data) as a Function of the Density-Height 
Sampling Interval, and as a Function of the Value of the Temperature, for Each 
of Nine Values of Relative Uncertainty in the Density-Height Observations 
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only a small influence. One may recall that these restrictions are as follows: 

1. The height h associated with T is more than 3 scale heights below 
hj the greatest height of the sounding, i.e. , (hj - h q ) > 3 H. 

2. The relative uncertainty of the density-height data is constant over 
the entire height range of the sounding. 

3. The atmosphere is isothermal over the height range of the sounding. 

The first of these restrictions must be retained unless one has some independent 
means for determining T, or Hj and its uncertainty as required by equation (33). 
Only for (h x - h q ) > 3 H will the ratio P 1 /p C{ be sufficiently small to permit the 
associated term to be neglected in equations (34) through (38). 

The second restriction is not a very significant one since it has been shown that 
the constant factor (6p/p), used as the general expression for density uncertainty 
in equation (35), is dominated by the particular value (6p q /p q ) associated with the 

f 

density-height data at h q . Thus, the need for (§p/p) to be constant over the 
entire height range of a particular sounding can be relaxed without significantly 
affecting the validity of any of equations (36), (37), or (38). 

The third restriction appears to be philosophically important because the earth’s 
atmosphere below 120km is certainly not isothermal. Actually, however, the 
exi stance of the various non- zero gradients in the temperature profile has little 
effect on 5T q /T q . This situation is due to the fact that ST q /T q depends pri- 
marily on the conditions between the data points (p q-1 , h q-1 ) and (p q , h q ), and 
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hardly at all upon the condition between other pairs of data points. Even a ndn- 
zero temperature gradient between the two specified data points has only a small 
effect, and this can be readily accounted for. This is accomplished by reintro- 
ducing into equation (36) certain temperature -gradient-dependent coefficients 
which are associated with Sp q /p q in the general version of 5H q /H q (or equiva- 
lently ST /T ), as expressed by equation (23). 

H M 

To begin with, (H/Ah) in equation (36) is replaced by its original form H/ (Ah + AH) 
as used in equation (23). Then, S q>min , which has been shown to be almost 
exactly equal to [(p - P q -\)/p q ] for isothermal conditions in the small-increment 
regime, is replaced by the product of this ratio times its associated coefficient 
a q _j from equation (23). Thus, for an atmosphere with varying temperature 
or scale height, we may rewrite equation (38) for the small -increment regime as 



Using the value of a as defined by equation (27), and replacing H by (2H q + 

AH)/2, where AH = (H ql - H q ), and remembering that Ah = (h q _ t - h q ), we 
have an expression for ST q /T q in terms of H q _ 1 and p q-1 at height h q-1 , and 
in terms of H q , P q , and 5p q at height h q : 



40 



It is interesting to examine an example of the use of equation (40) in terms of a 
realistic set of data, applicable to the small- increment regime, as taken from 
the U.S. Standard Atmosphere 1962. Choosing the height interval. Ah ~ 1km', 
between h q _ { = 66 km' and h q - 65 km', in a region of negative temperature grad- 
ient, AT/ Ah - -4 K/kmi we find p q-1 and p q to be 1.3482 x 10‘ 4 and 1.5331 x 
10" 4 kg m" 3 respectively, H q = 7.0709 km', and AH ~ -0.1173km'. The substi- 
tution of these data into equation (40) leads to a value of 1.076 for the coefficient 
of 6p q /p q , a value involving only a second-order difference from the unity coef- 
ficient of equation (38); i.e. , 5T q /T q , in this case, would be 1.076 times 6 P q /p q , 
instead of 1.000 times Sp a /p a for the isothermal case. Thus, the assumption of 
isothermality, both in the development and use of equations 36, 37, and 38, is 
seen to introduce only a second-order error, and accounting for this error is un- 
necessary in most uncertainty determinations. 
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V. CONCLUSIONS: A COMPARISON OF THE DENSITY-DATA NOMOGRAPH 
WITH THAT FOR PRESSURE DATA 

The nomograph relating 5T/T to sampling-height interval of density data for 
various temperatures and density uncertainties depicted in Figure 8 is based on 
the same range of values of temperature and sampling-height interval used in 
the nomograph of Figure 5, which relates 5T/T to the sampling -height in- 
terval of pressure data. In addition, the nine assumed values of 8p/p in Figure 
8 are identical to the nine assumed values of 6 p/p in Figure 5. Thus, the two 
figures are exactly comparable, with the only difference being that Figure 8 is 
for density data, while Figure 5 is for pressure data. A comparison of these 
two figures shows that, at least from the point of view of height resolution and 
uncertainty of derived temperatures, density -height data are to be preferred 
over pressure-height data. 

In the large -increment regime, i. e. , Ah > 5H, the nomograph of Figure 8 is 
essentially identical to that of Figure 5. In this regime, bothfigures show that 
the uncertainty in the derived temperature is decreased as the height resolution 
of the related temperature profile is decreased, i.e. , as the sampling-height 
interval is increased. Even the smaller sampling-height intervals wi thin this 
regime are so lahge, however, i.e.. Ah - 5H, or about 30km' , that the height 
resolution makes' this regime essentially useless for temperature -height de- 
terminations at heights below 120 km' . 
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In the small-increment regime as well as in the transition regime, Figure 5 and 


8 are quite different from each other. In the small-increment regime (Ah < 

0. 2H), the uncertainty in the temperatures derived from density -height data is 
seen to be independent of both Ah and T (or H), as long as (h - h ) > 3H. In 
this regime, 5T _/T derived from densities is dependent only upon the un- 

4 H 

certainty of the density -height value associated with the related height h q . It 
is apparent that uncertainty considerations place no limits on the usable fineness 
of the height resolution of the density -height data. This situation is in contrast 
to that associated with pres sure -height data in Figure 5 where ST/T is seen to 
depend upon T and Ah, as well as upon 5p/p, and where 5T/T becomes pro- 
hibitively large for reasonable values of <5 p/p when Ah becomes smaller than 
about 1km. 

In the transition regime, Figure 8 shows less difference from Figure 5 than in 
the small -increment regime. In this transition regime, the characteristics of 
the function determining 5T/T from density -height data vary between those of 
the small-increment regime and those of the large -increment regime, such 
that 6T/T decreases slightly as Ah varies from about 0. 2H to about 5H. In 
going from the small-increment regime to the transition regime, however, the 
increased coarseness of the height resolution of the related temperature -height 
profile would more than offset the correspondingly small decrease in tempera- 
ture uncertainty, particularly since the minimum values of Ah in this regime 
are already of the order of lkm'. Even in this regime, however, the 
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density -height data yield smaller values of 5T/T than are obtained from 
pres sure -height data. 

In general, a comparison of Figure 8 with Figure 5 shows that the density -height 
data are far more desirable than pressure -height data at least from the point 
of view of the size of the uncertainty and of the height resolution of the derived 
temperature -height profile. 
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APPENDIX A 


DEVELOPMENT OF THE GENERAL EXPRESSION FOR 
(5T q /Tq) 2 AND FOR (SH q /H q ) 2 


It is well known that temperature may be deduced from density -height data by 
means of the following relationship: 


P t GM 1 rh 

T = — • T + — - • — • / 

q Pq 1 * Pql 


P (h) dh 


(A-l) 


where 

P(h) is the function relating atmospheric density to geopotential height, 
hj is the geopotential height of the upper limit of the region of integration, 
in geopotential meters (m' ), 

hq is the geopotential height of the lower limit of the region of integration, 
in geopotential meters (m' ), 

G is the geopotential gravity constant 9. 80665 m 2 sec“ 2 (m' ) -1 , 

M is the mean molecular weight of the air 28. 9644 kg (kmol) -1 , 

R is the universal gas constant 8.31432 x 10 3 joules K _1 kmol* 1 , 

Pj is the atmospheric density at h., 

Tj is the atmospheric temperature at In, 

p is the atmospheric density at h , and 

Tqis the sought-after atmospheric temperature at h . 
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Because of the defined relationship between temperature and scale height H, 
i. e. , H = TR/GM, it is convenient to rewrite equation (A-l) as 


where ' • 

H 1 is the scale height at h x and 

H is the scale height at h . 

4 H. 


Pj 1 f h q 

H q “ — • H, + -•/ p (h) dh (A- 

4 P q p qJ h . ' 


Since p( h) is generally not known as an analytical function, for which one might 
find a perfect integral, but rather is known as a set of numerical values of 
density versus geometric or geopotential height, it is convenient to replace the 
integral term in equation (A-2) by a series approximation. One possible form 
of such an approximation is a series of terms each representing the area of 
successive trapezoids under the graph of the natural logarithm of density versus 
h. When the data, p versus h, is plotted on a semilogarithmic scale, i. e. , 
fin p versus h, with successive data points connected by straight-line segments, 
as is frequently the situation for closely spaced density -height data, the series 
of logarithmic trapezoids exactly fits the area under the graph. In using this 
approximation to the integral equation, (A- 2) may be rewritten as 


_ p, i A (Vi ytPj Pj-t) 

4 " Pq ' * Pq £ 


(A-3) 


The validity of the approximation of equation (A -2) by equation (A-3) improves as 
the height increment between successive density values decreases. 
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The uncertainty in the computed value of H is based on a function involving the 
partial derivative of H q with respect to each of the independent variables. These 
include H, and the appropriate number of density -height data pairs; i. e. ; h ) , 

Pj > hj , pj (for j=2toq-l);as well as h q , p q . With the assumption that the 
uncertainty in each data pair is entirely in the density value, we are interested 
in the partial derivatives of H q , as expressed by equation (A -3), with respect 
to only the following variables: H 1# p, , p, (for j = 2 to q - 1), and p , having 
the general designation y . . These partial derivatives dH / by- multiplied by 
the corresponding uncertainty dy i are given in the following equations; 


The product of times the partial derivative of H q with respect to H : is 
simply 


3H ( 

3H~ 


6 Hl 

1 Vq 1 J H . 


(A-4) 


The product of Bp 1 times the partial derivative of H q with respect to p ^ is 

!2a. 5p Jh.h ( ■ Ah Ufi ( J 

3p l 1 Vq 1 \£np 2 - Cn Pl 7 [_p q \,p q J \Stii Pl - tnpj \ J Pl (A ' 5) 


£np 2 - £np, 

It can be shown, however, that 

1 


Crip- - Strip 


- = (— ) 

! \ Ah + AH J j . 


(A-6) 


2 ""ri ■*•—**/, 2 

where the double subscript "1, 2" on the right-hand member of equation (A-6) 
indicates that Ah = (hj - h 2 ), AH = (H t - H 2 ), and H = (Hj + H 2 )/2. When 
is greater than h 2 » the left-hand side of equation (A— 6) is positive. Ah is 
positive, and Ah has the sign of the temperature gradient 3T/3h in the region 
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h to h ->♦ The same equation, with other consecutive digits, as for example, 
1 2 

2,3 or 3,4, etc. , indicates these same relationships for the corresponding 
height intervals. ' •- 


Equation (A-6) combined with equation (A -5) yields 


3H. 


— ~ * 8 p, 

a pj Pl 



Pl ( P 2-Pi\( H \ 11^1 

Pq \ Pq / \ Ah + AH /l,2jJ Pi 


(A-7) 


The product of 5 p 2 times the derivative of H q with respect to p 2 is a quantity 
which, when modified by the introduction of equation (A -6), with each of two 
different, but appropriate pairs of digits, becomes 



In this equation the subscript "1, 2 ,f on each of two factors has the same signi- 
ficance as it has on these same two factors in equation (A-7), while the subscript 
"2,3" on each of two other factors signifies that the quantities H, AH, and Ah 
within each of these factors are associated with the geopotential height incre- 
ment h 2 to h 3 . Thus, for these two factors, Ah “ h 2 - h 3 , H = (H 2 + H 3 )/2, and 
AH = H 2 -H 3 . Equation (A -8) obviously involves two height increments, the ,, 
two which are separated by the height h 2 . 


The partial derivative of H q with respect to each of p 3 , p 4 , . . . , p q _ 2 > and 
p , when multiplied by 5pi (where j is successively 3, 4, ...» q-2, and q-1) 

n 1 J 



and, when further modified by the appropriate introduction of equation (A -6), is 
identical to equation (A -8) except for the successive incrementing of the sub- 
scripts. Thus, a convenient form of the product of 5p q-1 , times the particular 
expression for the partial derivative of H q with respect to p ql is 



In equation (A -9) the subscript "q-2, q-l M on each of two factors signifies that 
the quantities H, AH, and Ah within each of those two factors are associated with 
the geopotential-height increment h q2 to h j , so that Ah = h 2 - h ql , H - 
(H q _ 2 + H q_ 1 )/2, and AH = H q _ 2 “ H q-i * Similarly, in that equation, the sub- 
script r 'q-l,q"on each of two other factors signifies that the quantities H, AH, 
and Ah within each of these two factors are associated with the geopotential -height 
increment h q l to h q , so that Ah = h q _j- h q , H = (Hq-x + H q )/2, and AH = H q-1 - 
H . Again this expression involves two height increments, in this instance the 

H. 

two separated by the height h , . 

4 1 

Because of the common form of the partial derivatives of H q with respect to 
p 2 , p 3 » . . . , p q _ 2 > and p , it is desirable to write a general version of equa- 
tions (A -8) and (A -9) to express the product of times the partial derivative 
of H with respect to p. , where j is understood to have values ranging from 
2 to q-1. 
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This general equation is 



Finally, the product of 5p q times the partial derivative of H q with respect to p q 
when modified by the introduction of equation (A -6), is 



In this equation the subscript on each of two factors has the same signi- 

ficance as in equation (A-9), while the subscript j" on another factor in the 
general term signifies that the quantities H, AH, and Ah within that tei’m are 
associated with the general geopotential height increment hj_j to hj . Thus, for 
this factor. Ah = hj_j - hj , H = (Hj_ 1 + Hj )/2, and AH = Hj_ x - H- . 

It is evident that one part of this equation deals with the data associated with the 
single height increment Ah = (h q _ x - hj between the lowest two density -height 
values involved in the calculation of H , while another part of the equation deals 

SI 

with all the height increments between h 1 and h^, and their associated density 
data. 
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Each of equations (A-4), (A-7), (A-10), and (A-ll) is directly involved in the 
Gausian expression for relative uncertainty 5H /H fl which follows: 

4 4 



The first, second, and fourth terms within the brackets of this equation are 
seen to represent exactly the squares of equation (A-4), (A -7), and (A-ll), 
respectively. The third term represents the squares of equation (A-10) evalu- 
ated for j ranging from 2 to q-1. The bracketed portion of the right-hand side 
of equation (A -12) is seen to be multiplied by the reciprocal of H , thereby im- 
plying that each term of equations (A-4), (A-7), (A-10), and (A-ll) must ulti- 
mately be multiplied by l/H q , either before or after these equations are squared 
and summed to equal (5H q /H q ) 2 . It is convenient in this case to do the multipli- 
cation before squaring and summing. Equations (A-7), (A-10), and (A-ll), each 
contain the doubly subscripted factor 

/H # Ah \ 

\ 1 * Ah + AH / 

in at least one term. It is convenient therefore to accomplish this multiplication 
operation in the appropriate terms by introducing H q into the denominator of this 
doubly subscripted factor, thereby converting this factor into a nondimensional 
coefficient with a value close to unity. Then, in order to put the resulting coef- 
ficient into proper prospective with respect to the remainder of the equation, 
and also to conserve space in the uncertainty expression being developed, it is 
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convenient to define each of these modified factors by a specific symbol. Thus, 
from equation (A -7) the modified factor is defined as the coefficient a 12 , i, e. , 

■ (A ' 13) 


while the two modified factors involved in equation (A -11) are defined as 

V.j = (h^ • Ah + AH^j j ' (A:14) 

and - •. 


yjj+i 


<x jj+ 1 ” f H q * Ah + AH ^ (A ‘ 15) 


The modified factor from equation (A-10) is defined as the coefficient q , 


i. e. , 


a 


q- 


= fjL. -_Ah_\ 

i,q |^H q Ah + AHy q l ( 


(A- 16) 


Substituting equations (A -4), (A -7), (A-10), and (A-ll), respectively, into the 
successive terms on the right-hand side of equation (A -12), dividing each of these 
equations by H q (i. e. , replacing each of the "l's" in the denominators of these 
equations by H q ), and simultaneously replacing the resulting modified doubly 
subscripted factors in these equations by the equivalent coefficient forms defined 
in equations (A-13) through (A -16), leads to the following expression for 


<«H q /H q )2: 
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Because 5H q /H q is identically equal to ST a /T a , equation (A-17) may be re- 


9 q 1 


written as 



J ,2 


P 


q \ 'q /V fih+A a/l.a 


2 \2 


{“ q "' q fe H 


- n / H 


p q j \ Ah + AH , 
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